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Abstract. The outbursts of low mass X-ray binaries are prolonged relative to those of dwarf nova cataclysmic variables as 
a consequence of X-ray irradiation of the disc. We show that the time-scale of the decay light curve and its luminosity at a 
characteristic time are linked to the radius of the accretion disc. Hence a good X-ray light curve permits two independent 
estimates of the disc radius. In the case of the millisecond pulsars SAX J1808.4-3658 and XTE J0929-314 the agreement 
between these estimates is very strong. Our analysis allows new determinations of distances and accretion disc radii. 
Our analysis will allow determination of accretion disc radii for sources in external galaxies, and hence constrain system 
parameters where other observational techniques are not possible. We also use the X-ray light curves to estimate the mass 
transfer rate. The broken exponential decay observed in the 2002 outburst of SAX J1808.4-3658 may be caused by the 
changing self-shadowing of the disc. 
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1. INTRODUCTION 

The outburst profile of X-ray accreting millisecond pul- 
sars (MSPs) shows a decay similar to those of tran- 
sient low mass X-ray binaries (LMXBs), i.e. exponen- 
tial or linear decays, with also a distinct knee. Only 
XTE J 1807-294 shows a purely exponential decay with 
a time scale of ^120 days. The outburst light curve of 
HETE J 1900. 1-2455 is completely different which have 
shown no exponential decays during the outburst @]. 
Compared to other accreting MSPs with outburst periods 
of a few days to a month, HETE J 1900. 1-2455 shows ev- 
idence of being a "quasi-persistent" X-ray source during 
the long outburst period [6], and references therein. 

King & Ritter (19989) (henceforth KR) and Shahbaz, 
Charles & King (1998) (henceforth SCK) examined the 
X-ray light curves of transient LMXBs in decay, char- 
acterising them as simple exponential or linear decays. 
In some cases, however, the decays show a distinct knee. 
As shown e.g. in Fig.|2]for SAX J 1808.4-3658 this is not 
a secondary maximum in the sense that the X-ray count 
rate declines throughout. 

In X-ray bright disc- accreting systems, including tran- 
sient LMXBs during outburst, the temperature of the ac- 
cretion disc is dominated by X-ray heating from the in- 
ner accretion regions across most of the disc. For a disc 
in which the scale height can be described as a function 
of radius by H = HqR" for constant Ho and n, the tem- 
perature, T, is given by lfl3ll as 

\%gR 2 R y ' 



where si£ is the disc albedo, Lx is the central X-ray lu- 
minosity heating the disc, assumed to be emanating ra- 
dially from a small spherical source. Only if n > 1 is the 
whole surface of the accretion disc illuminated; if n = 1 
Eq. (Q} predicts a surface temperature of zero which is 
invalid because of local viscous heat production. Consis- 
tent with KR we adopt a value of n = 9/7, although see 
section Modified exponential limit. Equation ([TJ can be 
rewritten to give the maximum radius which is heated 
by X-rays to the temperature 7), needed to remain in the 
hot, viscous state. Thus 

R " n =^ H « {n - l)L *- (2) 
We abbreviate this equation as 

^-"=4>Lx- 0) 

KR showed that X-ray heating during the decay from 
outburst causes the light curves of transient LMXBs to 
exhibit either exponential or linear declines depending 
on whether or not the luminosity is sufficient to keep the 
outer disc edge hot. They note that exponential decays 
must revert to the linear mode when the X-ray flux has 
decreased sufficiently, but do not analyse this in detail. 
Nor do they consider the effects of mass transfer from the 
donor star, — M2, during the outburst. Whilst in many sys- 
tems — M2 is negligible compared to the mass accretion 
rate onto the compact object, M c , during outburst, this 
is not necessarily the case. In this proceeding I review 
the results of the MSPs pubblished in Powell, Haswell & 
Falanga (2006). 



2. OUTBURSTS PROFILE 
2.1. Exponential decay 

The mass of the accretion disc during the exponential 
decay is given by equation (3) of KR; 



M di , 



3vkr 



(4) 



where Vkr is some measure of viscosity, taken by KR 
to be the viscosity near the outer disc edge with a value 
of Vkr ~ 10 11 m 2 s _1 . In the case of significant mass 
transfer into the disc from the donor star the central 
accretion rate is given by 



M c = -M 2 -M di! 



(5) 



Using this to replace M c in the mass of the hot disc 
can be written as 



Mdisc = M a exp 



(6) 



where M a is the constant of integration, and the X-ray 
luminosity is proportional to 



Mr 



-M 2 - 
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At some time, t t , the temperature of the outer disc edge, 
while still dominated by X-ray heating, will be only just 
sufficient to remain in the hot, viscous state, i.e. = 
^disc- We denote the corresponding X-ray luminosity 
as L t , with the central accretion rate M,. The X-ray 
luminosity at earlier times is 



Lx = (L,-L 2 )exp[ - 



3 Vkr (t-t,) 



R 2 

"disc 



L 2 , (8) 



where L 2 = 77 (— M 2 ) c 2 , and 77 ~ 0.2 is the efficiency 
with which rest mass is liberated as X-rays. 



2.2. Linear decay 

When /?/, < 7?(jisc the outer part of the disc is no longer 
kept in the hot, viscous state by central irradiation. The 
radius of the hot part of the disc will decrease with the 
decreasing X-ray luminosity, corresponding to the late 
linear decline of KR. This, too, may be modified by the 
donor star mass loss rate, depending on the mass ac- 
cretion rate through the now cold outer disc. Designat- 
ing fl c (R) as the rate at which mass is transported in- 
wards through the cold disc, we examine the two extreme 
behaviours of jj. c . For the minimally efficient cold disc 



jJ. c (R) = for all R, and none of — M 2 reaches the hot disc 
when Rh < ^disc- At the other extreme ju c -(/f) = — M 2 and 
the cold disc does not increase in surface density with 
time. The occurrence of outbursts suggests that the latter 
case is false; with stable cold-state mass transfer, out- 
bursts would never occur. 

The mass in the hot disc changes due to three factors: 
accretion from its inner edge, the loss of material into the 
encroaching cold disc, and mass transfer from the cold 
disc. Hence, 

M hot = -M c + Llc(Rh)+27CR h Z(R h )R h . (9) 

Using the approximate substitution from ([3]) that 

R\ - 4>T]c 2 M 6 ., (10) 

(O becomes 

M hot = -M c + n c {R h ) + K<Pric 2 -L(R h )M c . (11) 

Additionally, equation (2) of KR allows Mh ot to be writ- 
ten as 



M hot = 

Equating CD) to d/df ([HI 
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(12) 



1 



M C R\ = -M c + fl c (R h ) + n<t>T]c 2 Z(R h )M c , 



3vKRdf 
so using equation (2) of KR 
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M c 



(13) 



(14) 



We now substitute the two cases of — ju c (/?/, ) outlined 
above. When mass transfer into the hot disc is negligible, 



M c 



3 Vkr 

4>T7C 2 ' 



(15) 



and the decay is the linear decline predicted by KR as 
expected. When ll c (Rh) = —M 2 , 



M c 



3 Vkr (-M 2 -M c ) 
<Pric 2 M c 



(16) 



While M c 3> — M 2 this reduces to the same form. As 
M c decreases it will approach the limit of M c = —M 2 , 
at which point M c = 0. The exact form of the decay 
will depend on the unreliably known form of ^l c (R), and 
during the late linear decline viscous heating near R^ 
becomes non-negligible. We do not predict the form of 
the late decay other than to note that if there is some 
radius R\\ m satisfying 



M(*iim)T?C 2 = ^, 



(17) 



then the l.h.s. of (fTTT i will give the asymptotic limit of the 
decay and the accretion luminosity at the beginning of 
quiescence. 

As a final point on the linear decay we note that the 
viscous instability may cause a range of annuli within the 
cold disc to enter the hot state and rapidly transfer mass 
inward. Some of this may enter the inner hot disc, or else 
may increase the surface density of the subsequent cold 
disc such that the inner hot disc receives more material 
at its outer edge. Consequently, small rebrightenings or 
other artefacts potentially complicate the linear decay. 



2.3. Continuous derivative 

To examine the observed X-ray light curves of tran- 
sient LMXBs it is useful to examine the constraints on 
the transition. The gradient of the exponential decay at 
time t t is 

Lx = -^(L t -r 1 (-M 2 )c 2 ), (18) 

"disc 

whilst the gradient of the linear decline is given by ( fT31 ), 

3Vkr 



L X = — 
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From ^ we adopt 



(19) 



(20) 



where it is assumed that at the outer disc edge H 
0.27?disc- We therefore obtain 



Lx(tt) = - 



3vkr / n (-M 2 ) , 



U 



(21) 



In the case that — 77 M 2 c 2 is small relative to Lx the gra- 
dient of the exponential decay at time t t is equal to that 
of the subsequent linear decline; the first derivative of the 
X-ray light curve is expected to vary smoothly through- 
out the decay. If however, — 77 M 2 c 2 cannot be neglected 
then the linear decline is steeper than the exponential 
decay at t t and a discontinuity in gradient is expected. 
Whether this is detectable depends on the value of — M 2 
and the quality of the X-ray data. 

Some X-ray light curves, including most of those dis- 
cussed in this paper, include a knee feature in the decay 
from outburst. Where this can be explained by an expo- 
nential to linear transition with non-negligible — 77 M 2 c 2 
we label it as a 'brink' . In general the decay light curve 
may contain multiple knees due to different effects, but 
we do not expect more than one brink unless the source 
rebrightens between them such that each occurs at ap- 
proximately the same luminosity, corresponding to a disc 
not varying greatly in radius. 



3. APPLICATION TO OBSERVED 
X-RAY TRANSIENTS 

The full source samples which fits the outburst decay 
light curves predicted with the model described above 
can be found in Powell, Haswell & Falanga (2006). 
Here I report only the results of some of the MSPs 
sources. This work could be extended also to the MSPs 
IGR J00291+5934 H, XTE J1814-338 & and Swift 
J1756.9-2508 [16]. We use our analysis to determine 
distances and accretion disc radii. System parameters 
from the literature for the objects analysed are collected 
in Table Q] for comparison with quantities we derive. 
Since we will be calculating the radius of the accretion 
disc we calculate from the system parameters the orbital 
separation a. 

The count rates N as a function of time t were fitted 
with 

N = (N, -A^) exp f-^) +Ne (22) 

for the exponential phase, where N e is the limit of the 
exponential decay, is the count rate at the brink, and 
T e is the time-scale of the decay. The linear decline is 
given by 

N = N t (\- t -^j, (23) 

where T; is the time after the brink at which the count rate 
would become zero if the linear decline continued. Using 
the assumed distance to each source with a hydrogen 
column density of 10 22 cm~ 2 it is possible to convert the 
count rates Nt and A^ into absolute X-ray luminosities 
L t and L e , the X-ray luminosity at the brink and the 
limit of the exponential decay (corresponding to — M 2 ) 
respectively. For the systems, XTE J0929-314 for which 
we have simultaneous ASM and PCA light curves, the 
ratio of count rates is approximately 3 1 . We present the 
luminosities in the 1.5-12 keV range corresponding to 
the ASM spectral range using an assumed power law 
spectrum with photon index of T = 2. The luminosities 
thus determined are given in Table [2] We note that if 
the spectrum differs between systems our deduced L t is 
invalid, though our deduced exponential timescale will 
be unaffected. 

Using (0 and ([8]) it is possible to calculate the disc 
radius independently from the measured values of N t and 
T e ; the result based on x e is of particular importance as 
it does not depend on the conversion of count rate to 
flux but only on a directly observable timescale. It is 
also possible to determine — M 2 from N e . Results are 
given in Table for comparison we have also listed the 
circularization radius, ^ c iro an d the distance between the 
compact object and the inner Lagrange point, b\, for each 
source based on the values in Table Q] Discussion of the 




40 50 60 70 80 

Modified Julian Day - 52 400 



FIGURE 1. Top panel: The fitted ASM light curve of 
4U 1543-475. The solid line shows the exponential-then-linear 
fit, the dashed line indicates the best linear fit. The data shown 
in grey were not included in the fit, lying before and after the 
period of interest. Middle panel: Residuals to the linear fit. Bot- 
tom panel: Residuals to the exponential/linear fit. 

choice of values for 4> and Vkr is given in section Best- 
fitting parameters. 

3.1. An example: 4U 1543-475 

For compareson to the MSPs I show also the outburst 
profile from the black hole system 4U 1543-475 (see in 
Powell, Haswell & Falanga 2006 for more detail). Fig.Q] 
shows that the exponential-then-linear fit to the decay 
from the 2002 outburst of 4U 1543-475 is significantly 
better than the least-squares regression line. Because no 
discontinuity is seen in the gradient, the fitting routine 
used was based on the continuous derivative model (see 
section Continuous derivative). Previous outbursts were 
reported in 1971, 1983 and 1992. The source therefore 
spends only a small fraction of its time in outburst, sug- 
gesting that — M2 is small compared with the central ac- 
cretion rate during outburst, consistent with the source's 
failure to display a broken gradient. 

3.2. SAX J1808.4-3658 

SAX J1808.4-3658 was first detected in outburst in 
1996 and has subsequently shown outbursts in 1998, 
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FIGURE 2. The fitted PC A light curve of the 1998 outburst 
of XTE J1808.4-3658. The dashed line indicates the fit of SCK. 
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FIGURE 3. The fitted PCA light curve of the 2002 outburst 
of XTE J1808.4-3658 



2000, 2002, and 2005. Hence it is expected that the 
ratio of — Mi to the peak central accretion rate is much 
higher than in 4U 1543-475. This is consistent with the 
appearance of a broken decay in the two well observed 
outbursts, shown in Figs. [2] and [3] Each light curve is well 
fitted by a decay that is exponential towards a positive 
limit followed by a sharp gradient discontinuity (brink) 
leading to a linear decline. The fit to the 1998 outburst 
includes all points that can be confirmed to occur during 
the decay; the observation at day 16 of Fig.|2]may belong 
to a period of constant luminosity or even an unsteady 
rise to outburst, while the observations around day 33 of 
Fig. |2] may involve a rebrightening. Even if this is not 
the case the late linear decline is expected to become 
increasingly complex, as we noted in Section Linear 
decay. Comparison with Fig. [3] shows that the brink 
occurs at approximately the same count rate in both 
cases. This is expected if /?disc remains constant (i.e. L t 



TABLE 1. System parameters from the literature. Orbital separation a is based on a 1.4Mq neutron star. 



Source 


^orbit [min] 


Pspin [ms] 


a [10 8 m] 


M 2 [Af ] 


Distance [kpc] 


Ref. 


XTE J0929-314 


43.6 


5.4 


3.2 


0.008-0.03 


5 


\Si 


XTE J 175 1-305 


42 


2.3 


3.1 


0.013-0.035 




[17] 


XTE J 1807-294 


40 


5.25 


3.0 


0.01-0.022 




[2, 4] 


SAXJ1808.4-3658 


120.8 


2.5 


6.3 


0.04-0.1 


2.5 


[3, 29] 



TABLE 2. Decay light curve parameters. Luminosities are based on the photon index and distance 
from Table Q] and values of Nh from the same sources, and are stated for a spectral range of 1.5- 
12 keV. For consistency with Fig.|5]we take 1 ASM count s to equal 31 PCA counts s for our 
purposes, suggesting that T ~ 2.1. 



Decay 


L, [ lO^Js -1 ] 


L e [ KpJs -1 ] 


l e [ days ] 


%l [ days ] 


XTE J0929-314 


49±1 


48 ±1 


6.9±0.4 


19.7± 1.2 


XTE J1751-305 


361 ±1 


110±30 


5.9±0.5 


1.69 ±0.02 


SAXJ1808.4-3658 1998 


38.54±0.03 


32.3 ±0.2 


4.89 ±0.06 


2.978 ±0.007 


SAX J1808.4-3658 2002a 




50±7 


2.8±0.3 




SAX J1808.4-3658 2002b 


45.1 ±0.1 


29.6±0.6 


5.0±0.1 


3.22 ±0.03 



is the same). 

The dashed line in Fig. [2] indicates the fit proposed 
by SCK, in which they identified the surplus around the 
brink as a secondary maximum. We disagree and view 
their fit as inadequate in having a negative limit to the 
exponential decay. It is possible that features in other 
decay light curves identified as secondary maxima but 
having monotonically decreasing X-ray flux have the 
same explanation. 

A second issue in Fig. [3] is the broken nature of the 
light curve prior to the brink. Whilst we do not have a 
detailed explanation for this, we speculate that it may 
be associated with some regions of the disc being self- 
shadowed and initially remaining in the cool state. 3D 
hydrodynamic simulations e.g. [7] show self-shadowing 
is likely, due to spiral density waves and irradiated disc 
warping c.f. M25I1 . Consequently the effective area of the 
hot disc, the quantity probably corresponding to our mea- 
surement of R% sc , increases during the decay in an ap- 
parently stepwise manner, introducing new material into 
the hot disc and therefore resetting the central luminos- 
ity to a higher value. This view is supported by the fact 
that fitting the region around day 65 of Fig. [3] indicates 
a shorter exponential time-scale than from day 70 to the 
brink, corresponding to a smaller outer radius of the hot 
region. By day 70 we assume that the entire disc is in 
the hot state. The growing irradiated area could lead to 
an actual increase in Lx, hence it is an explanation for a 
type of secondary maximum. It is possible that the light 
curve will monotonically decline despite the growing ir- 
radiated area, hence producing knees in the light curve. 
The light curve of SAX J1808.4-3658 around day 67 of 
Fig.|3]is suggestive of this behaviour. 




FIGURE 4. The PCA light curve of the 2002 April outburst 
of XTE J1751-305, fitted after MJD 52 373. 

3.3. XTEJ1751-305 

XTE J 175 1-305 was first detected during its 2002 out- 
burst, shown in Fig. [4] and a second outburst was ob- 
served in 2005. If this repetition interval is normal for the 
system, it indicates a ratio of — Mi to peak central accre- 
tion rate more similar to that of SAX J1808.4-3658 than 
to 4U 1543-475. Consistent with this, the 2002 decay 
light curve has a shape similar to those of SAX J 1808.4- 
3658 in that a smooth exponential decay is broken by 
occasional rebrightenings. Therefore only the last por- 
tion, appearing to correspond to a single smooth decay, 
is fitted. The linear portion of this decay terminates at a 
constant level. In accordance with dT6b we interpret this 
as the rate at which mass can be transferred through the 
mostly cold disc; this rate will continue during quies- 



TABLE 3. System radii and mass transfer rate. The two values of -R(ji sc and —M2 are calculated from the X-ray 
light curves; i? c ; lc and b\ are calculated from the data in Table Q] Errors exclude the effect of the uncertainty in 
source distance. 



Decay 



"circ 

[ ioV; 



[ io 8 m; 



^disc(M 
[ 10 8 m] 



[ ioV; 



-Mi 
[ 10 1 n V/ . yr ] 



XTE J0929-314 


1.5- 


-1.9 


2.63- 


-2.82 


2.52 


±0.03 


2.67 ±0.08 


84±2 


XTE J 175 1-305 


1.4- 


-1.7 


2.5- 


-2.7 


6.85 


±0.01 


2.47 ±0.11 


200 ±50 


SAXJ1808.4-3658 1998 


2.06- 


-2.71 


4.68- 


-5.07 


2.238 


±0.001 


2.252 ±0.014 


57±4 


SAX J1808.4-3658 2002a 


2.06- 


-2.71 


4.68- 


-5.07 






1.70±0.09 


90 ±10 


SAX J1808.4-3658 2002b 


2.06- 


-2.71 


4.68- 


-5.07 


2.421 


±0.003 


2.28 ±0.02 


51.8±1.1 
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FIGURE 5. The fitted PCA light curve of XTE J0929-314. 
The ASM light curve, suitably scaled, is overplotted in grey. 



cence whilst the disc mass increases. 



3.4. XTEJ0929-314 

A third observed example of the brink is found in the 
only observed (2002) outburst of XTE J0929-314. Here 
the overall shape before the brink is exponential, but a 
maximum at about day 10 of Fig. [5] and a minimum at 
about day 19 are anomalous. Possibly the outer edge 
of the disc remained cold until late in the exponential 
decay and several secondary maxima of the type already 
discussed have occurred. 



3.5. XTEJ1807-294 

XTE J1807-294 is, like SAX J1808.4-3658, 
XTE J175 1-305 and XTE J0929-314, one of eight 
known millisecond pulsar (MSP) systems. MSPs have 
short orbital periods and hence small and therefore 
relatively simple discs. We therefore expect the X-ray 
light curve of XTE J 1 807-294 to also display the brinked 
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FIGURE 6. The PCA light curve of XTE J1807-294 

decay seen in the first three. Fig. [6] does not appear 
to lend itself to this type of fitting, and we have been 
unable to obtain reasonable fitting parameters. One 
possible explanation is that the linear decay had already 
begun at the start of the observed light curve; in this 
case the brink count rate must be around fifty per cent 
greater than that seen in SAX J1808.4-3658. The orbital 
separation of XTE J 1807-294 (c.f. TableQ} is half that of 
SAX J1808. 4-3658, so the disc must be correspondingly 
smaller. A high value of N t can only be accounted for 
if XTE J 1807-294 is at a distance of less than about 
1 kpc so the luminosity remains low while the countrate 
is high. Conversely, if the light curve corresponds to 
the exponential portion with several rebrightenings, the 
maximum value of may be 15, with the light curve 
after day 80 tentatively assumed to be linear. If the 
disc is half the radius of that of SAX J 1808.4-3658, 
in proportion to the orbital separation, (fTTT i indicates 
that the brink luminosity will be one quarter that of 
SAX J1808.4-3658. The PCA count rate at the proposed 
brink is one quarter that of SAX J1808.4-3658, indicat- 
ing a similar source distance of D ~ 2 kpc. We deem 
the latter explanation more likely, and suspect that the 
exponential decay of XTE J 1807-294 is complicated by 



self-shadowing of the disc. 

The more extreme mass ratio of XTE J 1807-294 than 
SAX J1808.4-3658, indicated in Table Q] suggests that 
the accretion disc radius in XTE J 1807-294 is a larger 
fraction of the orbital separation, indicating a larger brink 
luminosity than assumed above, increasing the distance 
estimates in both cases. A lower value of N t is possible 
in the second case, which would also indicate a greater 
source distance. 



4. BEST-FITTING PARAMETERS 

To force the disc radii based on T e to lie within the 
range R c i m < /?disc( T e) < b for the two best described 
systems, SAX J 1808.4-3658 and XTE J0929-314 the 
viscosity of Vkr = 4x 10 m 2 s was adopted, close 
to the value of 10 u m 2 s _1 proposed by KR. To make 
Rii sc {L t ) consistent with /?disc(Te) for most of the neutron 
star systems a value of <t> = 1.3 x 10 -12 m 2 sJ~ 1 was 
adopted. Since R^c is based on <PL t , our systematic 
underestimate of L t by considering only the flux in the 
1.5-12 keV band suggests a yet smaller value of 4>. 
Calculated disc radii using these parameter values are 
given in Table [3] 

The large value of R<a %c (L t ) for XTE J 175 1-305 may 
be the result of spectral differences, or in this case may 
be the result of an erroneous source distance estimate. 
XTE J 175 1-305 has been assumed to lie at a distance of 
8kpc,being possibly associated with the Galactic Centre 
[e.g. |9J]; if the true distance is more like 3kpc then the 
two values of Ra sc will be approximately equal. 

The increasing value of /?disc(f e ) for the 2002 decay 
of SAX J1808.4-3658 and the erratically increasing lu- 
minosity between days 65 and 70 of Fig. |3]are both con- 
sistent with the interpretation that the outer disc edge was 
not in the hot state at the beginning of the decay. As such, 
the disc radius measured corresponds at first to a hot disc 
smaller than the disc as a whole, which is therefore de- 
pleted on a shorter timescale. At some point the outer 
disc is heated and enters the hot state. As well as increas- 
ing the timescale of the decay, this introduces extra ma- 
terial into the hot disc, causing rebrightening. 



5. MODIFIED EXPONENTIAL LIMIT 

We have assumed in the initial model that the limit of 
the exponential decay is equal to the assumed constant 
mass transfer rate from the donor star. This requires that 
the hot disc is capable of supporting the additional mass 
flux — M2 at all times and radii. However, the exponential 
nature of the decay itself implies that as the hot disc is 
depleted in surface density, its mass transfer rate falls. 
The mass transfer equation of the hot disc is linear in 



surface density, which can therefore be decomposed into 
two components; one representing the decay analogous 
to an unfed disc, the other corresponding to the effects 
of feeding the hot disc at the rate —M2- The former 
component we denote L e and the latter £2. Because the 
viscosity of the accretion disc depends on temperature, 
and the temperature is declining during the decay, it is 
necessary for E2 to increase with time if it is to sustain 
the same mass flux, implying that the corresponding flux 
jj.2 must decrease inwards. Thus the limit of the central 
mass accretion rate does not correspond to — M%, but 
rather /^(^m). The relevance of this issue is made clear 
by the light curves of 4U 1705-44, in which the limit of 
the exponential decay is often below the long term mean 
count rate of 12 ASM counts s _1 (see (author?) JH] for 
details). 

The standard derivation of the diffusion equation [e.g. 
24] can be used to give the mass transfer rate inward 
through the disc jx as 

H=6nR l/2 (vLR l/2 y . (24) 
Using the relation 



li 1 = InKL, 



(25) 



we derive equation (9) of [15] as expected. We require 
the form given in (1241 to prevent a constant of integration 
arising. 

To make use of d24l i we must adopt an appropriate 
form for the viscosity v. We use the a-viscosity form 



v = ac s H, 



(26) 



where c s is the sound speed, and the scale height is given 
by equation (2) of lfl5ll : 



1/2 



(27) 



Using ([T} to substitute for T and assuming that H(R) is a 
power law we obtain 



H = 



\GM\m J V 47TC7 ' x 



for a mean particle mass m ~ m p . We write 
of constant Hq as 



H — H Q R dif 



U 



1/7 



R 



(28) 
in terms 



(29) 



where n = 9/7 is derived and Ho ~ 0.2 is adopted from 
KR. |l|] note that if the X-ray flux is assumed to be ab- 
sorbed at the mid-plane, the power law index is changed 



to 45/38; we therefore consider values of 45/38 < n < 
9 jl as valid to simplify subsequent derivation. Using this 
form for H in (Q3 to derive the viscosity we find that 



Lx\ 2/1 f_R_ 



v = v [ ^ | [ — I , (30) 



>-3)/8 



where Vo is defined as 



v = a 



AKOm A 



Ho (3D 



We simplify ( f30b by adopting n = 1 1/9, so ( 1241 becomes 



^disc \ Lt J 



If we impose pi = — M2 for all radii we obtain 

£2 



-M 2 f L x \ 2/7 fl d isc 



37TVQ V U 



R ' 



(32) 



(33) 



corresponding to a total contribution to mass in the hot 
disc of 



-2M 2 /LxV 2/7 ,2 



(34) 



where the radius of the inner disc edge is much smaller 
than /?disc- As Lx decreases, a greater surface density and 
total mass is required to transport the same — M2 inward, 
so some of this material must be absorbed to increase the 
surface density. For a typical luminosity at the beginning 
of the exponential phase 2.5 times greater than L t the 
change in mass required is 



A»J2 



2 ^disc 

3v 



(l-2.5~ 2/7 ). (35) 



Examining 4U 1705-44 (see (author?) JH for de- 
tails), for which we have the largest mean luminos- 
ity relative to L t , we find a mean long term count 
rate of 12 ASMcountss -1 , whereas the limit of the de- 
cay a for instance is 8.2ASMcountss~'. Taken over 
an exponential phase duration of 57 days this makes 
the Lh.s. of ([35j equivalent to 1.9 x 10 7 ASM units of 
mass. Using the conversion factor from WebPIMMS of 
1 ASM count s -1 = 4 x 10~ 13 Js~' m~ 2 , at a distance of 
7.3 kpc this integrated flux is equivalent to 3 x 10 20 kg. 
Using Vkr = 4 x 10 10 m _2 s _1 (as in Tabled we find 



^dis 



:6.4 x 10 8 m, 



(36) 



in strong agreement with the values given in Table [3] 
despite the weakness of the derivation. In principle it is 
necessary to simultaneously derive the time and radius 
dependencies of jj. given the outer boundary condition of 

M(fldisc) = -M2- 



6. CONCLUSIONS 

A knee in the light curve of the decay from outburst of an 
transient LMXB is a natural consequence of mass trans- 
fer onto the outer edge of the disc, since this supply is ef- 
fectively cut off from the compact object when the outer 
disc enters the cool low-viscosity state. When the knee 
can be interpreted in this way we refer to it as a brink, 
since knees of other types may be seen in the decay out- 
burst. The X-ray luminosity at which the brink occurs is 
that at which the outer disc edge is just kept hot by cen- 
tral illumination, allowing this radius to be calculated. In 
addition, the exponential time-scale of the decay gives a 
second measure of the disc radius; these two estimates 
are in good agreement. By examining systems with well 
constrained disc radii we deduced values of the constants 
<t> and Vkr required to perform the calculations. This al- 
lows more accurate results to be obtained in other sys- 
tems. Where the source distance is poorly constrained, 
the timescale x e allows the absolute luminosity L t to be 
estimated, providing a measure of source distance. 

As well as 'secondary maxima' associated with the 
brink, we have identified maxima whose natural expla- 
nation is the discontinuous increase in the radius of the 
hot disc during outburst. This increase arises from ini- 
tially self-shadowed disc regions becoming irradiated. 
This process can occur until the whole disc is in the hot 
state. It is characteristic for the time-scale of the expo- 
nential decay to increase with each such maximum, and 
it is expected that the limit of each exponential decay will 
also increase. Although we do not have sufficient data to 
examine the matter in detail, it is probable that some sec- 
ondary 'maxima' of this type will also involve monotonic 
decay, and therefore probably will appear as a knee in the 
X-ray light curve. 

Current instrumentation is capable of measuring 
the light curves of X-ray transients in M31, for 
example ll28"ll show part of a 2004 July decay of 
XMM J0043 15.5+412440. While this light curve, with 
only 4 data points over 3 consecutive days, is too sparse 
to allow us to deduce any parameters from fitting the 
decay, our method could be applied to better-sampled 
extragalactic transient decays. Since the distance is 
relatively well-known for objects in external galaxies, 
we could use our method to deduce their accretion disc 
radii. This would contribute to constraints on the fun- 
damental system parameters, such as orbital separation 
and mass ratio, which might otherwise be impossible to 
determine. 
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